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ABSTRACT 

Supermassive black holes are investigated as possible sources for low-frequency bursts 
of gravity waves. The event rate for 'known' supermassive black holes at intermediate and 
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high redshifts, inferred from the quasar luminosity function is low ~ 0.1 yr _1 . A number 
Q density of gravitational-wave sources comparable to the number density of galaxies inferred 

^ from faint galaxy counts is necessary to raise the event rate significantly above one per 
year. A space-based interferometer could therefore only see several events per year from 
Qh supermassive black holes if an additional population of supermassive black holes existed 
O and emitted gravitational waves efficiently. These might reside in the population of dwarf 

^ galaxies or in a transient population of small dark matter haloes that have mostly merged 

into larger haloes hosting the galaxies seen today, as proposed in a hierarchical cosmogony. 
In the latter case, event rates could be of order 10 — 1000 per year due to coalescing 
supermassive-black-hole binaries formed during the merging process. Event rates could be 
as high as a few per second if the dark matter in galaxy haloes consists of supermassive 
black holes in the mass range 10 3 — 10 6 Mq, and produced gravitational waves efficiently. 
The proposed space-based gravitational-wave interferometer LISA/SAGITTARIUS should 
detect most gravitational-wave events involving supermassive black holes above 10 4 Mq out 
to redshifts of z ~ 100. 

1 INTRODUCTION 

Bursts from supermassive black holes are — apart from periodic gravitational waves 
from known binary stars — the major possible source of low-frequency gravitational waves 
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in the frequency range 10 — 1Hz. A considerable amount of work has been done on 
determining the strength, polarization and waveform of gravitational-wave bursts and on 
extracting these signals from noisy data (Peters & Mathews 1963; Estabrook & Wahlquist 
1975; Armstrong, Estabrook & Wahlquist 1987, Thorne 1987; Tinto & Armstrong 1991 ; 
Schutz 1992; Thorne 1993), but little attention has been paid to the possible event rates 
since the first study of Thorne & Braginsky (1976). 

The existence of supermassive black holes is still not proven beyond doubt, but a pop- 
ulation of supermassive black holes in the mass range of 10 6 - 10 9 M Q is believed to be the 
only viable explanation for most of the observed activity in active galactic nuclei (Lynden- 
Bell 1969, Blandford & Rees 1992). Thorne & Braginsky found an event rate ~ O.lyr -1 
for a source density comparable to the space density of typical galaxies (as defined 
by a fit of a Schechter-function to the galaxy luminosity function) and speculated that 
the observed number density of 'active' galaxies indicating the existence of supermassive 
black holes might be considerable higher than that. Since Thorne & Braginsky 's paper 
the knowledge of the luminosity function of quasars and galaxies and of the build up of 
galactic-size structures has greatly improved. This renders a more detailed investigation 
possible. 

Quasar activity seems to be quite common at intermediate and high redshift and prob- 
ably most of the luminous galaxies have undergone a phase of quasar activity, but, despite 
this, the integrated number density of observed quasars is unlikely to exceed the number 
density of galaxies, even if the life time of the individual quasar is assumed to be rather 
short (Haehnelt & Rees 1993). A considerably higher number of supermassive black holes 
of more moderate mass (10 4 — 10 6 M©) could be contained in observed dwarf galaxies or in 
the building blocks of luminous galaxies which, in hierarchical cosmogonies, would be built 
up by merging of smaller units at moderate redshifts between z = 4 and z = 0.5 (White 
& Rees 1978). The large number of faint blue galaxies might be a direct observational 
indication for such a hierarchical build-up of structure involving high merging rates and 
probably a higher number of small galaxies than inferred from the local galaxy luminos- 
ity function (Tyson 1988, Guiderdoni & Rocca-Volmerange 1990, Kauffmann, White & 
Guiderdoni 1993). If the dynamically inferred dark matter in galactic haloes consists of a 
population of black holes, the fraction of all mass contained in black holes would be even 
a factor 10 6 higher than that inferred for the black hole population powering quasars. 

The second important question is whether a population of supermassive black holes is 
likely to emit detectable gravitational waves. 

In principle, a significant fraction of the rest mass energy could be released during the 
formation of the black hole if it involves a non-axisymmetric collapse or rotation phase, but 
the experience from calculations for stellar mass size sources shows that the efficiencies are 
likely to be small and in the case of continuous slow accretion they could be even negligible 
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(see Schutz 1989 for a review). 

One source of gravitational radiation would be a close binary containing two super- 
massive black holes. It seems possible that supermassive black holes are formed as binaries 
by fission of a bar (Rees 1984) but a more likely path for the formation of a supermassive- 
black-hole binary is via the merging of the nuclei of two dark-matter haloes each containing 
a black hole. The putative black holes would sink to the centre of the merged nuclei on a 
dynamical friction time scale forming an initially rather wide binary. The newly-formed 
nucleus would be gas-rich and accretion onto the binary could lead to a rapid reduction of 
the separation. Effective emission of gravitational radiation would finally lead to coales- 
cence (Begelman, Blandford & Rees 1980). 

Another possibility for the emission of gravitational waves by a supermassive black 
hole is a head-on collision with a compact stellar object or the capture of such an object 
into a relativistic orbit. 

Section 2 gives an analysis of which kind of gravitational-wave events involving su- 
permassive black holes could be seen by a space-based gravitational-wave interferometer 
using the specifications of two recent proposals (LISA/SAGITTARIUS) for such a laser 
interferometer discussed in section 3 (Danzmann et al. 1993, Hellings et al. 1993). Sec- 
tion 4 discusses the expected event rates for different scenarios and section 5 contains the 
conclusions. A Hubble constant of Ho = 50kms _1 Mpc -1 is used throughout the paper. 

2 SUPERMASSIVE BLACK HOLES AS SOURCES OF LOW- 
FREQUENCY GRAVITATIONAL WAVES 

The gravitational-wave amplitude produced by a burst releasing a fraction e of the 
rest-mass energy of an emitting object of total mass M is given by 



3 G(eMci) y/2 

ftburst- l 2?r2 c l R{z) 2 f J W 



(Thorne 1987), where 



t-)/ \ c l-q +q z-(l- qo)(l + 2q z) I 
H qfi(l + z) 

is the comoving coordinate distance as a function of redshift z. qo = 0/2 is the conventional 
deceleration parameter and equation (2) holds for a vanishing cosmological constant, A = 0. 
For an Einstein-de Sitter universe of critical density, O = pj p c ^ = 1, equation (1) becomes 
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Vrst = 1.1 x 10- 20 h^el'l (y^) (1 - (1 + ^)" 1/2 )- 1 , (3) 

where h$Q denotes the Hubble constant in units of 50 km s _1 Mpc -1 . The characteristic 
burst frequency is related to the mass of the black hole and the redshift at which the burst 
occurs as 



f c = [(3 3 / 2 GM/c 3 ) (1 + z)]- 1 = 3.9 (y^-) _1 (1 + z)- 1 Hz (4) 
(Press 1971). For a binary with component masses M\ and M 2} the amplitude is given by 



/ 9 \ 1/2 r* 5 / 3 ,, 

= 3.4 x 10- /,„ (^f^f ^ ,1 - (1 + ,)-/-',-. 



(Thorne 1987, Danzmann et al. 1993), where M = M x + M 2 and \i = MiM 2 /M are the 
total and the reduced mass and x = M 2 /M\ denotes the mass ratio. 

The detector sensitivity as a function of frequency determines the most and the least 
massive black hole (binary) for which the detection of a burst is possible. This depends 
of course on redshift, emission efficiency and, in case of a binary, on the mass ratio. The 
proposed space based laser interferometer LISA/SAGITARRIUS (see section 3 for details 
of the sensitivity limits) would see gravitational wave bursts in the wedge shaped area above 
and to the left of the thick solid curves in Fig. 1. This assumes e = 0.1 as might be realistic 
for the coalescence of a binary with equal masses, and the thin solid curves show the limits 
for e = 10 -3 and e = 10 -5 . For binaries of unequal masses the sensitivity is reduced by 
a factor x -1 / 2 at the frequency of coalescence and the thin solid curves can therefore also 
be interpreted as the limits for the coalescence of binaries with mass ratios x = 10 -2 , 
x = 10 -4 (this assumes that the form of the signal is known with sufficient accuracy to 
apply matched filtering). For the coalescence of binaries of small mass occurring at small 
redshift, the highest sensitivity is obtained long before coalescence due to the narrow range 
in frequency where the interferometer reaches its peak sensitivity. The thick dashed curve 
shows the lower mass limit for a 'wide' binary with equal masses and the thin dashed 
curves show the sensitivity for binaries with mass ratios x = 10 -2 , 10 -4 . Here the signal 
is assumed to be integrated for the last year before coalescence. 
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A gravitational- wave detector with the sensitivity of LISA/SAGITTARIUS would see 
most possible gravitational-wave events involving supermassive black holes above 10 4 Mq 
back to very high redshifts, thus maximizing the possible event rates. A detector with a 
somewhat lower sensitivity would still cover a considerable fraction of the (M, z) plane, 
but as will be shown in section 4, even a sampling out to high redshifts does not necessarily 
lead to event rates sufficient to guarantee a few events per year. 

Compact objects of 1 Mq in a relativistic orbit around a black hole of 10 6 Mq would 
be seen out to moderate redshifts z ~ 1 — 3. Events involving black holes of mass <J 10 3 
Mq are only seen at non-cosmological distances and a binary of stellar mass black holes 
(10 M©) at the distance of Virgo might be just visible. 

3 THE PROPERTIES OF LISA/SAGITTARIUS 

LISA (Laser Interferometer Space Antenna) is a proposed space-based Michelson laser 
interferometer in orbit around the Sun designed to detect low-frequency gravitational waves 
in the frequency range 10 -5 — 1 Hz not detectable on Earth due to unavoidable seismic 
noise (Danzmann et al. 1993). Due to the vastly increased armlength (5 X 10 6 km) com- 
pared to a ground based interferometer like LIGO, the comparatively moderate accuracy 
in the determination of changes in the arm length achievable in space is sufficient to reach 
a sensitivity similar to LIGO in its advanced form. The SAGITTARIUS project tries to 
reach the same sensitivity specifications with a space-based interferometer orbiting the 
Earth (Hellings et al. 1993). The armlength is similar (10 6 km) and the projected resid- 
ual spurious accelerations and position measuring noise are also comparable to the LISA 
project. 

The sensitivity will be limited by spurious accelerations at frequencies below 10 -3 Hz 
and by shot noise at higher frequencies. For randomly oriented sources a search for gravi- 
tational bursts in a year of data with a Wiener filter renders a sensitivity of 

^burst ~ WS h fl 12 (6) 

(Thorne 1987), where Sh is the spectral noise density of the detector (in Hz -1 / 2 ) and f c 
is the characteristic frequency of the burst. A corresponding one-year search for periodic 
sources gives 

Aperiodic ~ 10^(3.15 XlO^)- 1 / 2 , (7) 
5 



for a signal-to-noise ratio S/N = 5. The projected maximum spurious acceleration for 
LISA is 3 X 10 -13 (//10 -4 Hz) -1 / 3 cms -2 Hz -1 / 2 , and the planned limit for the accuracy of 
the armlength measurement due to frequency-independent shot noise is 2 X 10 -9 cm Hz -1 / 2 . 
For an armlength of 5 X 10 6 km, the sensitivity for bursts is then approximately given by 



1.3 x 10" 21 (//10- 3 Hz)" 11 / 6 , / < 10" 3 Hz, 

iburst ~ \ 1.3 x 10" 21 (//10- 3 Hz) 1 / 2 , 10" 3 Hz < / < 3 x 10" 2 Hz, (8) 

6.9 x 10" 21 (//3 x 10" 2 Hz) 3 / 2 , / > 3 x 10" 2 Hz, 



The break at 3 X 10 2 Hz is due to the fact that there A/2 equals the armlength of the 
interferometer. The corresponding sensitivity for periodic sources is 



7.1 x 10" 24 (//10- 3 Hz)" 7 / 3 , / < 10" 3 Hz, 

Aperiodic 

7.1 x 10" 24 (//3 x 10" 2 Hz,) />3xlO" 2 Hz. 



The sensitivity limits for SAGITTARIUS are almost identical. 

4 EVENT RATES FOR GRAVITATIONAL WAVES FROM SU- 
PERMASSIVE BLACK HOLES 

For a population of gravitational-wave sources with comoving density N s (z), emitting at 
a rate ^ em , the total rate of received gravitational-wave events originating around redshift 
z between comoving coordinate distance R(z) and R(z) + dR is given by 



dv rec (z) = v em (z) (1 + z)- 1 N B (z) 4vr R\z) ,^ w<p , f, (10) 

y/i - k (R(z)/ny 

where k = —1,0,1 for an open, flat or closed universe respectively, and 1Z = c/(Hq (1 — O — 
A/3Hq)). For an emission rate of one per time interval corresponding to a unit increment 
in redshift 

u em = Hq (1 + z) (O (1 + z) 3 + A./'SHq ) 1 ! 2 and with the transformation from comoving 
coordinate interval to redshift interval 

- k (R(z)/K) 2 = c H' 1 (O (1 + z f + A/3H 2 )- 1 / 2 dz equation (10) becomes 
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(11) 



For an Einstein-de Sitter universe this becomes 



= 5-6 ( Ns{z) A (!-(! + z)-W? yr~\ (12) 
dz 50 Vio-2Mpc- 3 / V V ; ; y y ' 

In the following I will assess the space density of possible sources and the corresponding 
event rates for an Einstein-de Sitter universe. 

4.1 Observed quasars as sources of gravitational waves 

Due to the large surveys for high redshift quasars undertaken in recent years, the quasar 
luminosity function and its evolution with redshift has been determined quite accurately 
out to a redshift z ~ 4 - 4.5 (Boyle et al. 1988; Boyle et al. 1991; Boyle 1993; Irwin, 
McMahon & Hazard 1991; Warren, Hewett & Osmer 1991; Warren, Hewett & Osmer 
1993). With the assumption that quasars are short-lived (life time £q ~ 10 8 yr) and that 
each of the black holes produces a gravitational wave event at the redshift where the quasar 
is observed (e.g. by coalescing of two black holes during, shortly before or shortly after 
the observed activity) the number density of observed quasars in a given luminosity range 
(L $(-L, z)) can be transformed into a corresponding number density of supermassive black 
holes emitting once per unit increment in redshift 



dA ^7» = ' , r „ L ML = SM). (13) 

dlnM t Q H (1 + zf/ 2 K ' K J 

$(L) dL is the number density of quasars in the luminosity range dL and a linear trans- 
formation of quasar luminosity to black hole mass 8 is assumed. The numerical value of 8 
is somewhat uncertain but is of secondary interest for an estimate of the event rates as a 
LISA-like detector would be sensitive to all events out to a redshift z ~ 100 for the mass 
range in question. Long- and short-lived models for the quasar population are discussed 
in the literature (e.g. Cavaliere et al. 1983), but a detailed examination of the quasar 
luminosity function and its evolution (Haehnelt & Rees 1993) seems to indicate a rather 
short lifetime of order of the Salpeter time tsalp ~ 4e ra( j ; o.l 10 7 yr (the e-folding time for 
the mass of a black hole accreting at the Eddington limit and radiating with efficiency 
e rad)- Fig. 2 shows the expected event rates for a given mass range as obtained from the 
quasar luminosity function integrated from redshift z = to redshift z\, where z\ varies 
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in steps of 0.5 (solid curves). The event rates would be dominated by quasars at redshifts 
z < 3 where the luminous function is rather unanimously determined and the luminosity 
function of Boyle et al. (1991) is used. At higher redshift the luminosity function is some- 
what controversial and in order to show that the assumed luminosity function at z > 3 has 
little effect on the total event rate, a constant space density is assumed for z > 3. This 
overestimates the number density of luminous quasars with z > 3 by a factor of up to 10 
compared to the one found by Warren et al. (1993). The differential rate (per redshift 
interval) is shown by the dashed curves. It is clearly seen that the number density of black 
holes inferred from the quasar luminosity function can not generate more than one event 
every few years if each black hole produces only one event during its active lifetime. Even 
a smooth extrapolation of the quasar luminosity function to smaller masses does not raise 
the rate above one per year. 

However, high redshift quasar surveys have a rather bright limiting absolute magni- 
tude and nothing is known about the faint end of the quasar luminosity function. An 
additional more numerous population of massive black holes of 10 4 — 10 6 Mq radiating at 
the Eddington limit would not have been detected in current surveys at redshifts larger 
than 0.5 — 1. It could furthermore be that such a population does not radiate efficiently. 

4.2 Faint/Dwarf galaxies as hosts of massive black holes 

In this section the event rates are investigated for the case that a population of massive 
black holes (10 4 — 10 6 Mq) resides in observed small galaxies and has produced gravitational 
waves efficiently at some redshift. To estimate the number density of sources, I first use 
the local galaxy luminosity function. Galaxy luminosity functions are well described by a 
Schechter function 



were a describes the slope at the faint end and are the characteristic number density 

and luminosity at the bend in the function. Fig. 3 shows the event rate inferred from the 
number density of all galaxies more luminous than a given luminosity L as a function of 
L/L* assuming 



for three different values of a (-1.07,-1.25,-2.0) and $* = 2 x /i^lO -3 Mpc" 3 . The thick 
solid curve (a = —1.07) corresponds to the luminosity function obtained by Efstathiou et 




(14) 




0.5 £ z £ zi + 0.5, 

otherwise, 



(15) 
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al. (1988) and a redshift for gravitational-wave emission z\ = 3.0. All galaxies down to 
0.01-L* are required to raise the event rate to a few per year. 

A considerably larger number density of small galaxies is inferred from deep CCD 
imaging (Tyson 1988). One interpretation of this apparent discrepancy between the faint 
galaxy counts and the faint end of the local galaxy luminosity function is a merging of 
these galaxies into larger galaxies (Guiderdoni & Rocca-Volmerange 1990, Guiderdoni & 
Rocca-Volmerange 1991, Kauffmann et al. 1994). With the assumption that each of these 
galaxies has emitted gravitational waves at redshift z\ (by formation of a supermassive- 
black-hole binary), the number counts of faint galaxies can be converted into a number 
density of gravitational-wave sources 



AT (xr \ ,47rAf/V(z ) z x - 0.5 £ z & z x + 0.5, 

W*o,*i,*)H Qtherwi (16) 



where M is the number of galaxies per steradian and V{zq ) the integrated comoving volume 
out to the typical redshift of the galaxies zq. Tyson et al (1988) find J\f = 4.5 X 10 5 deg -2 
for a limiting magnitude B = 27 and zq = 0.5. The obtained event rate is then somewhat 
more favourable 



and could be of order 30 per year (Fig. 4). 

4.3 Merging in a hierarchical cosmogony as source of gravitational waves 

The rather low possible event rate inferred for small galaxies is due to the flat nature of 
the local luminosity function at the faint end (for a = — 1 the number of objects per In L _1 
is constant). The faint counts already indicate that there might be more small structures, 
and it is a characteristic feature of hierarchical cosmogonies that they predict the number 

of haloes to rise roughly as M am+1 per InM -1 with a m ~ —1.8 1.9. For a constant 

mass-to-light ratio this would correspond to a luminosity function with a = a m close to 
—2 (see dotted curve in Fig. 3). In a hierarchical scenario, a specified dark matter halo of 
mass Mi has typically undergone a merging history involving ~ f3M\/Mo smaller haloes 
of mass Mo, where 1 — f3 indicates the fraction of mass accreted in the form of haloes less 
massive than M\ . Recent numerical and analytical studies of the merging of dark-matter 
haloes by Lacey & Cole (1993) indicate that 1 — f3 should rapidly decrease with increasing 
Mi /Mq, and I will consider f3 = 0.3 as a lower limit for big values of M\/Mq. If all building 
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blocks of mass Mo contained a black hole and emitted gravitational waves at redshift z\ 
the hierarchical build-up of galaxies brighter than would lead to an event rate 



V0.3/ VlO" 5 / V1O 4 M / y y } 

This assumes a mass of 5 X 10 11 Mq for an galaxy, a minimum black-hole mass M m [ n 
which can be seen out to very high redshifts, a ratio 7 e ff of black hole mass to mass of the 
dark-matter haloes merging and an emission at z\ > 1. 

However, the emission of gravitational- wave events due to the merging of putative black 
holes within merging dark-matter haloes does not necessarily have to be coupled to the 
occurrence of luminous galaxies. Hierarchical cosmogonies generally predict many more 
(small) dark-matter haloes than indicated by the number of (faint) galaxies (Kauffmann 
et al. 1994). For a given hierarchical cosmogony — characterized by its power spectrum 
of primeval density fluctuations — the Press-Schechter formalism allows simple analytical 
estimates of the number density of haloes as a function of mass and redshift 



dn halo 
dlnM 



2\V2 



7T 



El 

M 



<fln a 



dlnM 



6c 



(7 



exp 



(M,z) L l2a(M } z) 



-61 



(19) 



where cr(M, z) is the variance in the mass density po and 6 C ~ 1.68 is the linearly extrapo- 
lated overdensity at collapse (Press & Schechter 1974, Lacey & Cole 1993). Lacey & Cole 
give the merging rate of haloes of mass M with haloes of the same mass or larger 



'merg 



(M,z) /2\ 1/2 



H (l+zf/ 2 Vvr 



M 



a(M + M', z) (1 - a 2 (M + M> , z)/a 2 (M, z) fl 2 



(20) 



din a 



dln(M + M') 



exp 



A 2 



a 2 (M + M',z) a 2 (M,z) 



dM' 



M + M 1 



With the assumption that each merger leads to the formation of a supermassive-black-hole 
binary and subsequent emission of gravitational waves, this can be easily transformed into 
a number density of gravitational wave sources emitting once per unit increment in redshift 



dN s (M,z) _ dn halo ! 5/2 

dlnM " dlnM " mer§(M ' Z) ^ (1 + Z) ' (21) 



10 



The solid curves in figure 5 show the integrated event rates corresponding to mergers 
between redshift z and the present for two specific hierarchical cosmogonies, a standard 
cold dark matter scenario (with variance in the mass density in 8/ij~gg Mpc spheres as = 0.4) 
and a mixed dark matter model (25 per cent neutrino fraction, COBE normalization), but 
the numbers depend only weakly on the details of the chosen scenario. For halos smaller 
than 10 10 Mq the event rate exceeds 30 per year. If bursts from supermassive black holes 
with masses larger than 10 4 Mq could be seen out to redshifts z = 4, a black hole formation 
efficiency as small as = Mbi ac k hole/^halo ~ 10 -6 could guarantee a generous event 
rate for a space-based interferometer. Should the black hole formation efficiency be higher 
and should therefore the merging of smaller haloes lead to the emission of detectable 
gravitational-wave bursts, the event rates could be considerably higher. 

4.4 Black holes as dark matter in galaxy haloes 

One candidate for the dynamically inferred dark-matter haloes of galaxies is massive 
black holes in the mass range 10 3 - 10 6 M© (Carr, Bond k Arnett 1984, Hut k Rees 1993, 
Rix k Lake 1993). This would enhance the number density of black holes by a factor of 
order 10 3 — 10 6 compared to the estimates discussed in the previous sections. However, it 
is by no means clear that such a population of black holes would emit gravitational waves 
as the most likely mechanism — merging of the host objects and subsequent formation of 
a close binary — is missing. For the mass range in question, the event rate would only be 
sufficient if either the efficiency e is high enough to allow detection at high redshift or if 
the whole population is synchronized to undergo a gravitational-wave burst at very low 
redshift (Bond k Carr 1984). In the first case the event rate is given by 



where fi^h 1S the fraction of the critical density in black holes of mass M and emission at 
z\ > 1 is assumed. A small fraction of close binaries coalescing in a Hubble time would 
therefore be sufficient to reach an event rate above a few per year. 

4.5 Influence of the cosmological parameter and the redshift of the gravita- 
tional wave bursts 

For a given comoving number density N(z) of sources emitting once around redshift z\, the 
total received rate increases monotonically with growing z\ due to the increasing comoving 
volume sampled. The event rates increase by a factor 4.1 between z\ = 1 and z\ = 5, and a 
further increase does not change the event rates considerably. For z\ = oo the rate would 
be a factor 2.9 higher than at z\ =5. Lowering O and/or introducing a cosmological 
constant increases the comoving volume at a given redshift. For O = 0.2, A = 0, the 




(22) 
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event rate would be higher by a factor (1.4,3.2,25.0) at z = (l,5,oo) respectively and 
for O = 0.2, A = A/3Hq = 0.8, the corresponding factors are (2.0,3.1,3.8). However it 
should be kept in mind that the bigger comoving volume also reduces the estimates of the 
comoving density of gravitational-wave sources inferred from high redshift objects. 

5 CONCLUSIONS 

While a detector with the sensitivity of LISA should certainly see gravitational radiation 
emitted from binary stars in our or nearby galaxies, there are several assumptions required 
in the case of gravitational wave bursts from supermassive black holes: 

• Supermassive black holes have to exist. 

• Supermassive black holes have to emit gravitational waves effectively either when they 
are formed, or as (coalescing) supermassive binaries. 

• The number of gravitational-wave events involving supermassive black holes must be 
considerably higher than the one inferred from the number density of observed quasars. 
The number of sources has to approach or exceed the number of all (including the 
faintest) observed galaxies to guarantee event rates above one per year. 

The first assumption is widely accepted to be true and the second and the third would 
be very plausible in a hierarchical cosmogony if black hole formation occurs in the large 
number of small dark matter halos that inevitably have to form and a considerably fraction 
of which has to merge into larger structures. Rather inefficient black-hole formation in these 
small haloes would be sufficient, and typical event rates are of order 10 — 1000 yr _1 . 

The question of whether the formation of such a population of low-mass black holes is 
likely is difficult to assess. While it is certainly plausible that the gas in these haloes is 
able to cool and loose angular momentum (Loeb 1992), the potential wells would be rather 
shallow so that feedback processes of massive stars might prevent a concentration of the 
gas. On the other hand, the formation of one supermassive star of 10 4 Mq collapsing into 
a black hole would be sufficient. 

If the dynamically inferred dark matter haloes of galaxies consisted of black holes of 10 3 — 
10 6 Mq — a considerably more speculative assumption — , a rather small fraction of these 
holes producing gravitational waves efficiently (e.g. in binaries coalescing within a Hubble 
time) would guarantee a sufficient event rate. If all these black holes produced gravitational 
waves efficiently, the event rate could reach a few per second exceeding the confusion limit 
and producing a continuous background. 
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In summary, it can be said, that the population of supermassive black holes that a space- 
based gravitational-wave interferometer could see, has not yet been observed by other 
means. However, the existence of such a population is certainly plausible, and the fact 
that some speculation is involved might be outweighed by the fact that such a population 
could probably not be detected otherwise. 
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FIGURE CAPTIONS 



Fig. 1. The LISA sensitivity limits for gravitational waves emitted from supermassive black hole 
(binaries) of total mass M as a function of redshift. LISA is able to detect bursts with emission efficiency e 
in the the wedge shaped area to the left and above the solid curves. The dashed curves show the reduced 
lower mass limit for coalescing binaries with mass ratio x. The dotted curve shows the detection limits for 
the spiral-in of a 1 Mq compact object. 

Fig. 2. The rate of gravitational-wave bursts produced by black holes of mass M as inferred from 
the quasar luminosity function. For z < 3 the luminosity function determined by Boyle et al. (1991) is used 
with the assumption that each quasar produces one gravitational-wave event during a lifetime of 4 x 10 7 yr. 
For 3.0 < z < 4.5 the luminosity function is assumed to be constant. Solid curves show the integrated 
event rate between z = and z = (0.5,1.0,1.5.2.0,2.5,3.0,3.5,4.0,4.5) respectively. Dashed curves show 
the differential event rate per corresponding redshift interval. Event rates scale inversely with the assumed 
quasar lifetime. 

Fig. 3. The event rate inferred for galaxies more luminous than a certain fraction of from 
a Schechter-type galaxy luminosity function. It is assumed that each galaxy hosts a black hole capable of 
producing one detectable gravitational- wave event at redshift z. Solid curves are for the galaxy luminosity 
function as determined by Efstathiou et al. (1988) with a = -1.07 and $* = 2 x 10~ 3 /if Mpc -3 . The 
assumed emission redshift is z = 3 (thick curves) and z = 1.0,5.0 (thin curves). The dashed curve is for 
a = —1.25 and the dotted curve for a = —2.0. 

Fig. 4. The event rate for galaxies brighter than B = 27 inferred from faint counts (Tyson et al. 
1988) with the assumption that the average galaxy in the faint counts has redshift zo and each galaxy hosts 
a black hole capable of producing a detectable gravitational- wave event at redshift z\ . The event rate scales 
linearly with the assumed normalization (Af = 4.5 x 10 5 deg -2 ). 

Fig. 5. The integrated event rates for dark-matter haloes of mass M merging between redshift z and 
the present for a standard cold-dark-matter scenario (ag = 0.4) and a mixed-dark-matter model (neutrino 
fraction 25 per cent, COBE normalization). It is assumed that one gravitational- wave burst is emitted at 
the merging redshift. 
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